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Abstract. We report the first comprehensive investigation of the line profile variation caused by non-radial pulsation in a 
magnetic oscillating chemically peculiar star. Spectrum variation of the well-known roAp star HR3831 is detected using very 
high-resolution high signal-to-noise spectroscopic time-series observations and are followed through the whole rotation cycle 
of the star. In contrast to previous photometric pulsational measurements of roAp stars, our spectroscopic observational ma- 
terial admits straightforward astrophysical interpretation and, hence, opens new exciting possibilities for direct and accurate 
analysis of the roAp pulsations. We confirm outstanding diversity of pulsational behaviour of different lines in the HR3831 
spectrum and attribute this phenomenon to an interplay between extreme vertical chemical inhomogeneity of the HR3831 at- 
mosphere and a running pulsation wave, propagating towards the upper photospheric layers with increasing amplitude. The 
rapid profile variation of the Ndm 6145 A line, which shows the maximum pulsational disturbances in the studied wavelength 
region, is characterized by measuring changes of its equivalent width and the first three moments. Each of these observables 
exhibit almost purely sinusoidal variation through the pulsation cycle, with the amplitude and phase clearly modulated by the 
stellar rotation. We demonstrate that rotational modulation of the radial velocity oscillations cannot be fully explained by an 
oblique axisymmetric dipole (f = 1, m = 0) mode, implied by the classical oblique pulsator model of roAp stars. Pulsational 
variation of the higher moments, in particular the line width, reveal substantial contribution of the high-order it = 3) spherical 
harmonics which appear due to distortion of pulsations in HR 3831 by the global magnetic field. We interpret observations with 
a novel numerical model of the pulsational variation and rotational modulation of the line profile moments in roAp stars. The 
comparison between observed and computed amplitudes and phases of the radial velocity and line width variation is used to 
establish parameters of the oblique pulsator model of HR 383 1 . Furthermore, definite detection of pulsational variation in lines 
of light and iron-peak elements enables the first 3-D mapping of pulsations in non-radially oscillating star. 
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1 . Introduction ters of roAp pulsators (e.g., Matthews et al. U999l Cunha et al. 

1.1. RoAp stars Since ^ discovery of roAp pu i sat i ons by Kurtz JT9781 it 

became clear that strong magnetic fields in these stars have a 
defining role in exciting the oscillations and shaping the main 
pulsation properties. It was found that the amplitude and phase 
of the rapid light variation are modulated with the stellar ro- 
tation and that the phases of the magnetic field and pulsation 



The rapidly oscillating Ap (roAp) stars are cool magnetic 
chemically peculiar stars, pulsating in high-overtone non-radial 
acoustic p-modes. There are 34 roAp stars known at present 
time (Kurtz & Martinez 2000 for a recent review). These ob- 
jects oscillate with periods in the range of 6-21 min (Kurtz & ,. , „ ., , 
w i _„ . . i.,,,,-! amplitude extrema typically coincide with each other (Kurtz 
Martinez 2000, Elkin et al. 2005), while their light variation ,, A ,f., „ ., , , ^ 3 , „ , . , . 



1982 1. Detailed frequency analyses of the roAp photometric 
light curves convincingly demonstrated that main characteris- 
tics of their rapid brightness variation can be attributed to the 
oblique axisymmetric dipole (I = 1, m = 0) modes, aligned 
with the axis of nearly axisymmetric quasi-dipolar magnetic 

....... . . , field. This oblique pulsator model, first proposed by Kurtz 

Photometric investigations of roAp stars were carried out dur- , . , , . , . , , 

. ... (1982), gave rather successful and straightforward phenomeno- 



amplitudes rarely exceed 10 mmag in Johnson B. Pulsations 
in roAp stars bear certain resemblance to the solar acoustic 
p-mode oscillations, but are characterized by smaller number 
of excited modes and have three order of magnitude larger 
amplitudes compared to the p-mode oscillations in the Sun. 



ing the last 25 years and have yielded unique asteroseismic in- 
formation on the internal structure and fundamental parame 



logical explanation of the main features in the roAp frequency 
spectra. However, subsequent studies of roAp pulsations re- 

vealed that the mode geometry in some stars defied a simple in- 

* Based on observations obtained at the European Southern terpr etation in terms of a single spherical harmonic (e.g., Kurtz 
Observatory, La Silla, Chile (ESO program No. 66.D-0241 A). et al. 1989 1. Moreover, gradual accumulation of the photomet- 
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ric information about pulsations in some best-studied monope- 
riodic roAp stars with oblique dipole modes provided evidence 
for significant deviations of the mode geometry from the ex- 
pected axisymmetric dipole shape - an effect loosely called 
"distorted dipole modes" (Kurtz et al. 1997 1. 

Several theoretical investigations (e.g., Dziembowski & 
Goode [T9851 Shibahashi & Takata IT9"9"3l Cunha & Gough 
2000 1 have argued that combined effects of the magnetic field 
and rotation may be responsible for the departure of the sur- 
face geometry of pulsational perturbations from a pure spheri- 
cal harmonic. In fact, one of the most recent theoretical studies, 
presented by Bigot & Dziembowski ( 2002 1, has challenged the 
general framework of the classical oblique pulsator model of 
Kurtz ( 1982 1 by showing that accurate non-perturbative treat- 
ment of the interaction between stellar rotation, pulsations, and 
magnetic field may lead to an entirely new geometrical picture 
of the roAp oscillations. Calculations by Bigot & Dziembowski 
(2002 1 suggested that pulsation eigenmodes are represented 
by a complex superposition of spherical harmonic functions 
and may contain substantial non-axisymmetric components. 
Furthermore, strong effects due to the centrifugal force are ex- 
pected to break the locking of the pulsation axis with the dipo- 
lar magnetic field, so that pulsations are aligned with neither 
the magnetic nor rotation axes. It should be noted, however, 
that calculations of Bigot & Dziembowski ( 2002 1 are limited 
to the p-mode pulsations in stars with magnetic field strength 
^ 1 kG and are not directly applicable to the majority of roAp 
pulsators, which have surface magnetic fields of a few kG. 

In contrast to the approach followed by Bigot & 
Dziembowski J2002> . Saio & Gautschy ( I20041 and Saio 
( 2005 1 studied the effects of strong magnetic field on axisym- 
metric p-mode oscillations assuming an alignment between 
pulsations and dipolar magnetic field and disregarding distor- 
tion of pulsation modes by the stellar rotation. Saio & Gautschy 
( 2004 1 find substantial deviation of the surface structure of pul- 
sational fluctuations from a single spherical harmonic. In par- 
ticular, the I — 1 mode is distorted by a dipolar magnetic field 
in such a way that pulsations are described by a superposition 
of the spherical harmonics with odd angular degree I and pul- 
sation amplitude is strongly confined to the magnetic axis in 
the outer stellar layers. 

In the view of the diversity and increasing sophistication 
of the theoretical speculations about the geometry and physics 
of roAp oscillations it becomes apparent that modelling of the 
photometric light curves of roAp stars is generally insufficient 
for accurate determination of the structure of pulsation modes 
and is unable to provide useful tests of theoretical predictions. 
Information content of the high-speed photometric observa- 
tions is small due to averaging of pulsational disturbances over 
the visible stellar hemisphere (see Saio & Gautschy 2004 1 and 
highly uncertain because rapid light variation in roAp stars in- 
volves non-linear and non-adiabatic effects which are not at all 
understood (Watson IT9881 M edupe 120021 Saio 120051 It has 
been argued (e.g., Kurtz 1990 1 that spectroscopic time-resolved 
observations, in particular investigations of pulsational varia- 
tion in the Doppler-broadened spectral line profiles at differ- 
ent rotation phases, can potentially provide much more precise 
and unprecedentedly complete information about the vertical 



and horizontal structure of p-modes and its relation to the mag- 
netic field topology, chemical inhomogeneities, and anomalous 
atmospheric structure of Ap stars. Such spectroscopic studies 
will eventually allow us to address the underlying questions 
about the physics of pulsating peculiar stars, for instance, the 
problem of excitation of the magnetoacoustic oscillations or 
nature of the mode selection mechanism. 

1.2. Spectroscopic studies of the roAp pulsations 

Until recently our knowledge of the p-modes oscillations in 
roAp stars was based essentially on the high-speed photo- 
metric observations. Early time-resolved spectroscopic studies 
(Matthews et al. I1988I Libbrecht ll988> were focused on the 
detection of pulsational radial velocity (RV) changes and de- 
riving the ratio of RV and photometric amplitude. Subsequent 
more detailed low spectral resolution investigations of ffCir 
and HR3831 by Baldry et al. ( 1998bl [T998al revealed unex- 
pected diversity of the pulsational amplitude and phase in dif- 
ferent wavelength regions, showing that for roAp stars the con- 
cept of a "typical" amplitude of RV oscillations is meaningless. 
Baldry et al. ( 1999) and Baldry & Bedding (I5UU0T ) also found 
systematic change of pulsational amplitude and phase in the 
wings of the hydrogen Ha line, which was suggested to origi- 
nate from a strong depth-dependence of the pulsational charac- 
teristics in roAp atmospheres. 

High-resolution observations of y Equ by Kanaan & Hatzes 
(1998) localized RV variations of individual metal lines and 
confirmed that pulsations change dramatically from one spec- 
tral line to another. This phenomenon remained a mystery how- 
ever, because no clear underlying systematic trends or depen- 
dencies were identified. A breakthrough discovery was made 
by Savanov et al. ( 119991 ). who demonstrated that in yEqu the 
largest pulsational RV shifts of up to 1 kms -1 occur in the 
lines of Pr m and Nd m. Follow up very high spectral resolu- 
tion time-resolved observations of this star by Kochukhov & 
Ryabchikova ( 2001a) confirmed these results and also revealed 
the phase shifts between variation of singly and doubly ionized 
lines of the rare-earth elements (REE). Our study emphasized 
a stunning discrepancy (sometimes up to two orders of magni- 
tude in RV amplitude!) between strong pulsations in REE lines 
and the absence of changes in the lines of light and iron-peak 
elements, which do not show oscillations above few tens of 
ms _1 . Discovery of the pulsational spectroscopic variability in 
10 Aql (Kochukhov et al. 2002) and first analyses of the varia- 
tion of individual metal lines in a Cir and HR 3831 (Kochukhov 
& Ryabchikova 2001b) showed that this unusual pulsational 
behaviour, dominated by the REE lines, is not limited to y Equ, 
but is present in a very similar form in other roAp stars. This 
finding was further strengthened by the wide wavelength region 
analy ses of HR 1217 (Balona & Zim a l2002l . HR3831 (Balona 
I20021 ). a Cir (Balona & Lanev l2003l . 33 Lib (Mkrtichian et al. 
120031 and HD 166473 (Kurtz et al. l2003l . 

A plausible explanation of the nature of pulsational spectro- 
scopic variation of roAp stars was suggested by Ryabchikova 
et al. (2002). This study of the atmospheric structure of yEqu 
found compelling evidence for the strong vertical separation 
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(stratification) of chemical elements. Under the influence of 
radiative diffusion light and iron-peak elements are accumu- 
lated at the bottom of the y Equ atmosphere, whereas REEs are 
pushed high above normal line-forming region. Ryabchikova 
et al. ( 2002 ) proposed that the atmospheric layers contributing 
to the absorption in REE lines are located substantially higher 
compared to the formation depths of other elements. The most 
recent detailed non-LTE analysis by Mashonkina et al. ([2005 ) 
and Ryabchikova ( 2005 1 indicates that the Nd n and Nd m lines 
in yEqu originate from the optical depths 10~ 4 -10~ 5 . These 
high atmospheric layers rich in REE are characterized by en- 
hanced amplitude of non-radial pulsations. Within the frame- 
work of this model, the phase shifts between variation of dif- 
ferent lines are attributed to a running magneto-acoustic wave 
propagating outwards in the stellar atmosphere and increasing 
in amplitude with height. Similarities in spectroscopic pulsa- 
tional behaviour of different roAp stars suggest that they share 
similar vertical chemical and pulsational profiles. The method 
of tomographic mapping of the vertical structure of roAp pul- 
sation modes, pioneered by Ryabchikova et al. (2002i and 
Kochukhov ( 2003 ), uses vertical chemical inhomogeneities as 
spatial depth filters and aims at reconstructing detailed vertical 
cross-section of pulsational perturbations in roAp atmospheres, 
as can be done for no other type of main sequence pulsating 
star. 

Despite a dramatic recent progress in our understanding of 
the vertical structure of the roAp pulsation modes and the inter- 
play between pulsations and vertical stratification of chemical 
elements, relatively little attention has been paid to the prob- 
lem of inferring the horizontal geometry of pulsations in roAp 
stars. The question of systematic mode identification, central 
to the studies of other types of pulsating stars, has not been 
thoroughly investigated in the case of magnetic pulsators. To 
a large extent such an unsatisfactory situation is explained by 
the absence of suitable observational material. Studies of the 
vertical structure of pulsating cavity capitalize on the extreme 
vertical inhomogeneity of the roAp atmospheres and typically 
make use of a relatively low S /N spectroscopic time series, 
sufficient just to achieve a detection of pulsations and merely 
deduce amplitudes and phases of the RV changes in individ- 
ual spectral lines. In contrast, probing the horizontal structure 
of pulsation modes is necessarily based on the observations of 
the pulsational line profile variation (LPV) and, therefore, re- 
quires spectroscopic data of outstanding quality. Such observa- 
tions have been presented so far only for one ro Ap star, y Equ 
(Kochukhov & Ryabchikova 2001ai. In the study of this ob- 
ject we investigated variability of the Prm 6160 A and Ndm 
6145 A lines and probed the horizontal structure of the mag- 
netoacoustic /?-mode oscillation using several mode identifica- 
tion methods. Most importantly, it was convincingly shown that 
the moment variation of the REE lines is inconsistent with a 
picture expected for axisymmetric (to = 0) modes of low de- 
gree. This discovery of non-axisymmetric pulsations in a roAp 
star gave certain credibility to the generalized oblique pulsator 
model of Bigot & Dziembowski ( 2002 1 and demonstrated enor- 
mous diagnostic potential of the high-precision time-resolved 
spectroscopy of roAp stars. 



Despite its interesting pulsational properties, y Equ is not 
an ideal target for mapping the horizontal mode structure 
because its extremely slow rotation does not allow to put 
useful constraints on the surface chemical inhomogeneities 
and explore oscillations from different geometrical aspects. 
Additional complication is introduced by multiperiodicity and 
short growth and decay time of the p-modes in y Equ (Martinez 
et al. fT996l . Evidently, a more fruitful investigation of the LPV 
due to non-radial oscillation can be carried out for a monope- 
riodic roAp star with a high amplitude of pulsational spectral 
variation and a reasonably short rotation period. The present 
paper concentrates on such an analysis of the well-known roAp 
star HR 3831 and presents the first comprehensive study of the 
pulsational LPV in a roAp star. 

1.3. HR3831 

HR3831 (HD 83368, IM Vel) is one of the best studied rapidly 
oscillating Ap pulsators. This star is a brighter component in a 
visual binary with a separation of 3.29" and a AV = 2.84 mag 
difference in brightness. HR3831 has a visual magnitude of 
V = 6.25 (Hurly & Warner |T$H3} and belongs to the SrCrEu 
group of Ap stars (Houck 1978 1. The mean longitudinal mag- 
netic field of HR 383 1 was measured by Thompson ( 1983 1 and 
Mathys ( I1991> . The longitudinal field varies sinusoidally with 
an amplitude of 737 G and a mean value close to zero. The rota- 
tion period of P rot = 2.85 1962 d was derived from the magnetic 
measurements. 

Spectroscopic studies of HR3831 were conducted by 
Polosukhina et al. ( 2000 1 and Shavrina et al. ( 1200 11 1. They anal- 
ysed abundances of several chemical elements and studied in 
detail strong rotational modulation of the profile of the Li i res- 
onance doublet. 

Kochukhov et al. (2003 2004) carried out the first abun- 
dance Doppler imaging of HR 3831 using average spectra from 
the present study and additional high-resolution observations. 
We combined Hipparcos parallax, Stromgren and HyS photo- 
metric indices, and results of the model atmosphere analysis 
of the hydrogen lines in HR 3831 to derive fundamental stellar 
parameters: T eS = 7650 K, log g = 4.2, log L/L Q = 1.091, 
R/R Q = 2.003, and inclination angle of the stellar rotation axis 
i = 68.2°. It was found that many lines in the HR3831 spec- 
trum show noticeable variation with rotation phase, implying 
that strong horizontal abundance gradients are present at the 
stellar surface. In particular, the Doppler inversions showed Li 
and Eu to be concentrated in pairs of opposite small spots, lo- 
cated at the stellar rotational equator, whereas O exhibits an 
extended overabundance area along the big circle perpendic- 
ular to the rotational equator. In the context of the present 
study, behaviour of the doubly ionized Nd and Pr lines is of 
special interest because these REE lines turn out to be ideal 
pulsation indicators and are extensively analysed in the time- 
resolved spectra of roAp stars. Investigation by Kochukhov et 
al. (|2004 1 demonstrated that Nd m and Pr m lines do not show 
strong profile anomalies except for the rotation phases close to 
magnetic crossover (jp = 0.25,0.75), when pronounced dou- 
bling is evident in both the Nd and Pr absorption features (see 
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also Kochukhov & Ryabchikova 2001b I. Our Doppler mapping 
of HR3831 attributed this phenomenon to a ring of relative 
REE underabundance - a surface distribution just the opposite 
to that of oxygen. 

The horizontal abundance structures observed in HR3831 
are understood to be a result of the radiatively driven diffu- 
sive separation of chemical species in the presence of a dipolar 
magnetic field inclined at a large angle B with respect to the 
rotation axis of the star. This magnetic topology is in perfect 
agreement with the parameters of the dipolar model (B = 86.8°, 
Bp = 2.49 kG) determined by Kochukhov et al. J2004L but con- 
tradicts the very large 1 1 kG mean quadratic field measured in 
HR3831 by Mathys ( 1995 1. However, the latter measurements 
are known to be severely biased by the neglect of surface chem- 
ical inhomogeneities in the quadratic field diagnostic procedure 
(see Sect. 4 in Kochukhov et al. l2004> . 

Non-radial pulsations in HR3831 were discovered and in- 
vestigated by Kurtz ( 1982 ». Subsequent analysis of the pulsa- 
tional behaviour using photometric data spanning four observ- 
ing seasons was presented by Kurtz et al. (Q997 1. They revealed 
strong rotational modulation of pulsations and determined a 
principal pulsation frequency v = 1428.0091 fiHz (1 1.67 min). 
The power spectrum around the main frequency is described by 
a septuplet structure with a separation of the frequency com- 
ponents by v ro t = 4058.265 nHz, corresponding to the stellar 
rotational period P rot = 2.851976 + 0.00003 d. Pulsational 
light variation of HR3831 is clearly non-sinusoidal with the 
first-harmonic quintuplet and the second-harmonic triplet un- 
ambigously detectable in the Fourier transform of the time- 
resolved photometric data. 

A preliminary discovery report of the pulsational varia- 
tions in individual metal spectral lines in HR3831 was pub- 
lished by Kochukhov & Ryabchikova (2001b I. In this study we 
demonstrated that, similar to the pulsations in other roAp stars, 
the doubly ionized REE lines in HR3831 are distinguished 
by high-amplitude RV variation during pulsation cycle of the 
star. This was confirmed by Balona (2002 1, who obtained high- 
resolution time-series observations in a wide spectral region at 
several rotation phases of HR 3831. Balona found that all metal 
lines with large pulsational variability are invariably associated 
with the Prin and Ndm spectral features. He also detected a 
strong decrease of the amplitude of the pulsational RV varia- 
tion from Hff to the higher member of the Balmer series, which 
supported the idea that the amplitude of non-radial oscillation 
is a strong function of height in the atmospheres of roAp stars. 

The quality of the observational material collected by 
Balona (2002 1 was not sufficient to detect and study pulsational 
changes in the profiles of individual spectral lines. Instead 
Balona attempted to use co-added cross-correlation function 
derived from the Nd m lines for tentative mode identification. 
He suggested that the two frequencies detected in the RV vari- 
ability can be attributed to the I — 2, m — and I — 2, m — -2 
non-radial modes aligned with the stellar rotation axis. This 
mode identification remains controversial since it was not sup- 
ported by any quantitative analysis (e.g., line profile calcula- 
tions) and clearly contradicts the well-known dipole oblique 
pulsator model required to explain rapid photometric variabil- 
ity of HR3831. 



Table 1. Observing log of the time-resolved spectroscopic 
monitoring of HR3831. Columns give the UT date, heliocen- 
tric Julian Date for the beginning and end of observation, the 
duration of continuous monitoring in hours, the range of rota- 
tion phases, calculated according to the ephemeris of Kurtz et 
al. (119971 . and the number of spectra obtained on each night. 



Date 


HJD-2451900 


AT 




N 


05/02/2001 


45.5533-45.8779 


7.79 


0.965-0.079 


300 


06/02/2001 


46.5568-46.8796 


7.75 


0.317-0.430 


300 


07/02/2001 


47.5470^17.8678 


7.70 


0.664-0.776 


300 


13/02/2001 


53.5328-53.8786 


8.30 


0.763-0.884 


320 


14/02/2001 


54.5286-54.8716 


8.23 


0.112-0.232 


320 


15/02/2001 


55.5209-55.8678 


8.33 


0.460-0.581 


320 



A new pulsation Doppler imaging method of resolving 
horizontal structure of stellar non-radial oscillations was de- 
veloped by Kochukhov (2004a). This technique uses high- 
resolution observations of LPV and permits direct mapping of 
the surface pulsation velocity amplitude without a priori pa- 
rameterization of pulsation with spherical harmonic functions. 
Preliminary results of the application of this promising method 
to HR3831 were reported by Kochukhov ( 2004b I. Doppler 
mapping of pulsations in HR3831 was based on the spectra 
analysed here and suggested that distorted dipole oscillation 
mode is inclined at a large angle with respect to the stellar rota- 
tion axis and is aligned with the dipole component of the stellar 
magnetic field. 

2. Observations and spectra reduction 

Very high-resolution spectroscopic time-series observations of 
HR3831 presented in this paper were collected in February 
2001 with the Very Long Camera of the Coude Echelle 
Spectrograph (CES) fibre-linked to the Cassegrain focus of 
the ESO 3.6-m telescope. The spectrograph was configured to 
reach a resolving power of A/AA = 123 500 with the help of 
the medium resolution CES image sheer. The detector, the ESO 
CCD-61 2Kx4K chip, allowed us to record a 40.4 A spectral 
interval centred at A 6136 A. 

During CES observations we continuously monitored 
HR3831 for about 8 hours on each of the 6 observing nights 
with a break of 5 nights in the middle of the run. Such an ob- 
serving strategy enabled us to obtain 300-320 stellar spectra 
each night and ensure a good coverage of the rotation period 
of HR 3831, with phase gaps not exceeding 8% of the rotation 
cycle (see a log of observations in Table A total of 1860 
spectra was collected for HR 3831. Observations were obtained 
with a 70 s exposure time, whereas approximately 20 s was spent 
on reading the CCD detector and file transfer. The dead time 
was minimized by exposing only half of the detector in the 
cross-dispersion direction and additionally binning by a factor 
of four in the same direction. This setup allowed us to observe 
HR3831 with a 78% duty cycle and has proved to be an ap- 
propriate compromise between a reasonable sampling of the 
pulsation period and reaching a sufficient signal-to-noise ratio 
in individual spectra. 60 s exposures of the ThAr comparison 
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Fig. 1. Comparison between the mean spectrum of HR 3831 (thin line) and magnetic spectrum synthesis calculations (thick line) 
in the 6116-6156 A wavelength region is illustrated in the upper panel. Identification is given for the strongest spectral lines. 
Length of the lines underlining identification is proportional to the strength of the corresponding spectral features. The lower 
panel shows the standard deviation spectrum. 



spectrum were obtained every 1.3—1.4 hours during observa- 
tions. 

A dedicated set of IDL computer routines was developed 
to perform optimal reduction of the CES spectra based on the 
strategy and procedures implemented in the reduce package 
by Piskunov & Valenti (2002i. The reduction steps included 
construction of the master flat field and bias frames, subtrac- 
tion of bias and correction of the target spectra by the nor- 
malized flat field, subtraction of background, extraction of ID 
spectra, and correction for the overall curvature of the spec- 
tra in the dispersion direction. Since the CES instrument was 
designed to record observations at resolving power of up to 
A/AA = 235 000, its instrumental profile is oversampled in the 
medium resolution mode (t»5 pixels per resolution element). 
We took advantage of this situation and resampled all spectra 
in the dispersion direction, descreasing the total number of pix- 
els by a factor of two. This procedure yielded ID spectra with 
S IN in the range of 1 10-150 per pixel. 

A special procedure was developed to establish the wave- 
length scale and to monitor possible instrumental drifts of the 
spectrograph during observing nights. We identified and mea- 
sured positions of 23 emission lines in all available ThAr 
spectra and then fitted the pixel-wavelength relation with a 
third-order polynomial. Sequencies of comparison spectra ob- 
tained on the same nights were analysed simultaneously to de- 
rive common higher order polynomial coefficients and a time- 
dependent instrumental shift (the lowest order polynomial co- 
efficient). The final internal accuracy of the wavelength scale 
established through this procedure was » 8 x 10~ 5 (4 m s ). We 
found that the zero point of the CES spectrograph exhibited an 



overall smooth drift of 60-180 ms~' due to changes of ambi- 
ent temperature and pressure during observations of HR 3831. 
At the same time, the average shift between the consecutive 
ThAr spectra taken 1.3-1.4 hours apart was at the level of 17- 
35 ms _1 . This confirms an excellent short-term stability of the 
CES instrument and suggests that possible random instrumen- 
tal shifts do not exceed ~ 10 ms _1 . Consequently, it was justi- 
fied to linearly interpolate the zero point between the compar- 
ison spectra taken before and after groups of stellar exposures. 
Finally, the spectra were shifted to the heliocentric reference 
frame. 

In the last reduction stage the individual spectra of HR 3831 
were normalized to the continuum. This procedure was carried 
out in a series of steps in order to ensure a consistent normal- 
ization of all stellar spectra. First, a third order polynomial was 
fitted through the continuum points of the nightly mean spectra. 
This average continuum profile was applied to individual ob- 
servations and was subsequently corrected for each spectrum 
by the heavily smoothed ratio of the average and individual 
spectra. A similar procedure was followed to achieve a homo- 
geneity of the continuum normalization of the spectra obtained 
on different nights. 

3. Spectrum synthesis and line identification 

The line identification in the 6116-6156 A wavelength region 
of HR3831 was carried out based on magnetic spectrum syn- 
thesis calculations. We adopted the model atmosphere param- 
eters r e ff = 7650 K, log g = 4.2 and the projected rotational 
velocity v e sin« = 33kms _1 determined by Kochukhov et al. 
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Fig.2. Illustration of rapid variation of the Ndm 6145.07 A line in the spectrum of HR3831. The upper panels show LPV 
recorded on the night of 5 Feb 2001 and the lower panels show observations obtained on 7 Feb 2001. Each of the six horizontal 
subpanels presents the average spectrum of HR 3831 at the specified rotation phase and time sequence of the difference between 
the average and 50 individual time-resolved observations obtained during 74 min and covering roughly 6 pulsation periods of 
HR3831. Consecutive difference spectra are shifted in the vertical direction. The bottom part of the horizontal subpanels shows 
the standard deviation for each pixel of the observed spectra. 



(20041. The relevant parameters of atomic transitions were ex- 
tracted from the VALD database (Kupka et al. 1999 1 and were 
further fine-tuned by fitting the solar atlas. Synthetic spec- 
tra were calculated using the polarized radiative transfer code 
synthmag (Piskunov 1999 1. A homogeneous magnetic distri- 
bution was adopted and a field strength of B s = 1 .75 kG was 
estimated from the dipolar magnetic field model derived for 
HR 383 1 by Kochukhov et al. ( 120041 . 

The comparison between the average spectrum of HR 3831 
and our theoretical calculation is illustrated in Fig. [2 In gen- 
eral, we find a satisfactory agreement. Remaining discrepan- 
cies are attributed to averaging over spotted surface abundace 



structures (not accounted for in our spectrum synthesis) and 
unidentified spectral lines missing in the theoretical spectrum. 

In the previous time-resolved observations of the 
A 6150 A region in slowly rotating roAp stars (Kochukhov & 
Ryabchikova 2001a, Kochukhov et al. 20021 we were able to 
study several dozen spectral lines. Unfortunately, a moderately 
rapid rotation of HR 3831 makes a detailed pulsation analysis 
possible only for a few relatively unblended spectral features 
and blends with well-defined components. Those included the 
Cai 6122.22, Cm 6129.22 A lines, the blend of Sii and Cen 
at A 6131.9 A, the Ban 6141.71, Ndm 6145.07, Fen 6147.74, 
and 6149.26 A lines. Note that the latter Fen line is severely 
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Fig.3. Illustration of the Ndm 6145.07 A line profile variation at different rotation phases of HR3831. Each of the greyscale 
images is based on groups of 23 time-resolved spectra covering roughly 3 oscillation cycles of HR3831 at specified rotation 
phases. The upper plot in each panel shows residuals from the mean line profile. The greyscale plots show time evolution of 
the residuals, phased with the main pulsation frequency v = 1427.9958 //Hz and using the scale of +2.5% of the continuum 
intensity. The six columns of panels show LPV recorded during six different observing nights. The UT date is specified at the 
top of each column. Evolution of the variability pattern during observing nights is illustrated by three subpanels in each column 
which display spectra obtained at the beginning, middle, and end of the continuous monitoring. 



blended with an unidentified REE line at A 6148.85 A which of yEqu (Kochukhov & Ryabchikova 2001ai and other sharp 
shows a high-amplitude pulsation variability in the spectrum line roAp stars. 
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4. Pulsational behaviour of the Nd in 6145.07 A line 

4.1. Observations of the rapid line profile variability 

The high quality of the spectroscopic observational data ob- 
tained for HR 3831 enables a direct detection and comprehen- 
sive analysis of the LPV in REE spectral lines. The strong line 
of doubly ionized neodymium at A 6145.07 A was found to 
exhibit rapid line profile changes which are clearly modulated 
with the stellar rotation. Fig. |2 illustrates variability of this 
Ndm line at two representative rotation phases. At the phase 
of magnetic extremum, tp = 0.02 (the upper panels in Fig. |2}, 
the average profile of the Nd feature is not significantly dis- 
torted by the inhomogeneous surface distribution of this ele- 
ment. Pulsational variability has a simple symmetric pattern, 
expected for the oblique dipole pulsation mode. The line pro- 
file as a whole becomes skewed to the red and then to the blue. 
This is very different from the pulsational spectroscopic be- 
haviour of some of the hot non-radial pulsators which is char- 
acterized by high-contrast bumps and dips moving across stel- 
lar line profiles (e.g. Schrijvers et al. 2004). Such a behaviour 
is not observed in HR 3831 because of the low angular degree 
I of the excited mode. The dominant pulsational contribution 
comes from the global I - 1-3 mode which distorts significant 
part of the line profile, not unlike the spectral changes typical 
of the radial stellar pulsation (e.g. Telting et al. 19971. 

In contrast to the phase of magnetic extremum, pulsational 
changes of the Nd m line at the crossover phase (lower pan- 
els in Fig. |3 have lower amplitude and are more complex. 
According to the surface Doppler maps of HR3831 published 
by Kochukhov et al. (2004), rotation phase ip = 0.72 corre- 
sponds to a moment when the area of relative Nd underabun- 
dance passes through the centre of the stellar disk. This ex- 
plains doubling of the average line profile and suggests that 
complexity of the oscillation behaviour at this phase is due to 
a transition from one pulsation pole to another. Due to a mod- 
erately rapid rotation of the star and the respective substantial 
Doppler broadening, the areas associated with both pulsation 
poles are resolved in the wavelength domain and contribute 
to variability in the blue (respectively red) components of the 
Ndm line. It is worthwhile to emphasize that, despite a low 
amplitude and complexity of the oscillation pattern observed at 
the crossover phases, pulsation does not vanish entirely and can 
still be readily detected with our high-precision observations. 

An alternative representation of the Nd m line variability is 
shown in Fig. [2 These greyscale plots encompassing the whole 
dataset illustrate pulsational changes of the residuals with re- 
spect to the mean spectrum for the beginning, middle, and end 
of each of the six observing nights. The maximum deviation 
from the average spectrum reaches +2.5% of the stellar con- 
tinuum. Pulsation pattern changes dramatically from one night 
to another. The data also show unambiguous evidence for the 
systematic evolution of the variability pattern during individ- 
ual observing nights, capturing rotational change of the aspect 
angle which becomes noticeable during a few hours of continu- 
ous monitoring (for instance, note an increase of variability on 
6 Feb 2001 and decrease on 14 Feb 2001). 



The spectrum variability plots presented in Figs. and 
establish a smooth pattern of the rapid line profile changes in 
HR3831, similar to the pulsational variability observed in the 
atmospheres of many other types of stars pulsating with longer 
periods and low amplitudes. This observation represents a non- 
trivial result in the light of the chemical stratification analysis 
of Ryabchikova et al. l|2p02 1. They have demonstrated that the 
strong Nd m lines populating the roAp spectra form in tenuous 
high atmospheric layers, whose physical conditions are rather 
uncertain and where some complex non-linear pulsation effects 
may in principle take place. However, our observations do not 
provide any evidence for violent phenomena, such as super- 
sonic shocks or large non-linearity, associated with the prop- 
agation of pulsation waves. Consequently, we believe that, in 
the context of the present study, the standard analysis of the 
pulsational spectroscopic observables, carried out under (of- 
ten implicit) assumption that line profile changes are domi- 
nated by linear pulsation velocity fluctuations is fully justified. 
Nevertheless, it is useful to keep in mind that the spectroscopic 
diagnostic explored in our paper is sensitive to the pulsational 
perturbations in the upper photosphere and the respective pul- 
sation geometry is not necessarily the same as for the broad- 
band luminosity fluctuations originating in the lower part of 
the roAp atmospheres. 

The data presented here is the first observation of the pulsa- 
tional LPV in a rapidly rotating roAp star. The variability pat- 
tern found for the Nd m 6145.07 A line at the rotation phase of 
maximum pulsation amplitude is close to the one expected for 
an oblique dipolar oscillation (Kochukhov 2004a). The wave 
in the residual spectra propagates symmetrically from the blue 
to red line wing and then backwards. This behaviour differs 
from profile changes of the same Ndm line in the spectrum 
of very slowly rotating roAp star yEqu. For the latter ob- 
ject Kochukhov & Ryabchikova (2001a) discovered asymmet- 
ric variability characterized by the blue-to-red moving features 
only. However, observations of HR 3831 at the rotation phases 
of the magnetic crossover respectively minimum pulsation am- 
plitude reveal a LPV pattern which is more similar to the pul- 
sational changes observed for y Equ. Variability is limited to a 
narrow wavelength region in one or both wings and the wave 
in the residual spectra propagates asymmetrically from blue to 
red (Fig. |3] 14 Feb) or from red to blue (7 Feb). 



4.2. RV measurements and frequency analysis 

A low angular degree of the non-radial mode excited in 
HR3831 is reflected in a fairly simple pulsational line pro- 
file variation. Therefore, the integral observables, like the low- 
order line profile moments (Aerts et al. 1992), are sufficient 
for the purpose of frequency analysis and should yield essen- 
tially the same information as pixel-by-pixel line profile mea- 
surements (e.g. Telting & Schrijvers 1997). We determined mo- 
ments of the Ndm 6145.07 A line using expressions given by 
Kochukhov & Ryabchikova (2001a) whose formulas are equiv- 
alent to the standard line profile moment definitions of Aerts et 
al.(f!992l. 
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Fig. 4. Radial velocity variation of the Ndm 6145.07 A line in the spectrum of HR 3831. Symbols show individual measurements 
phased with the rotation period P = 2.851976 d. The smooth curve illustrates a Fourier fit of the rotational modulation of the 
average RV. 



Table 2. Results of non-linear least-squares fit of the Nd m 6145.07 A line profile moment variations. 



Frequency 


Wi (M) 


(V) (kms- 1 ) 


(V 2 ) (kirrV 2 ) 


(V 3 ) (km 3 s -3 ) 




A <p 


A ip 


A ip 


A <p 



The fundamental septuplet 



V 




3v rol 


227 


± 


80 


0.242 


± 0.056 


0.125 


±0.013 


0.180 ±0.018 


11.07 


±0.39 


0.392 


± 


0.006 


120 


± 


18 


0.557 


± 


0.024 


V 




2v rol 


175 


± 


75 


0.625 


± 0.069 


0.012 ±0.016 


0.806 ±0.166 


0.75 


±0.37 


0.665 


± 


0.079 


56 


± 


17 


0.554 


± 


0.049 


V 




Vrot 


2643 


± 


74 


0.926 


± 0.004 


1.079 


±0.013 


0.054 ± 0.005 


8.30 


±0.36 


0.312 


± 


0.007 


877 


± 


17 


0.033 


± 


0.003 


V 






249 


± 


79 


0.420 


+ 0.051 


0.025 


±0.015 


0.292 ± 0.084 


0.73 


±0.39 


0.029 


± 


0.084 


56 


± 


18 


0.174 


± 


0.052 


V 


+ 


Vmt 


1887 


± 


74 


0.851 


± 0.006 


0.594 


±0.014 


0.983 ± 0.006 


22.21 


±0.36 


0.723 


± 


0.003 


175 


± 


17 


0.967 


± 


0.016 


V 


+ 


2v ro , 


265 


± 


75 


0.232 


± 0.045 


0.018 


±0.013 


0.550 ±0.165 


0.93 


±0.37 


0.293 


± 


0.063 


45 


± 


17 


0.805 


± 


0.061 


V 


+ 


3v ro , 


34 


± 


80 


0.132 


± 0.378 


0.051 


±0.013 


0.187 ±0.040 


3.81 


±0.39 


0.804 


± 


0.016 


141 


± 


18 


0.420 


± 


0.020 




















The 


first-harmonic quintuplet 






















2v 




2v rot 


118 


± 


74 


0.986 


±0.100 


0.038 


±0.012 


0.030 ±0.051 


1.34 


±0.36 


0.272 


± 


0.043 


3 


± 


17 


0.441 


± 


0.775 


2v 




Vrot 


66 


± 


74 


0.808 


±0.178 


0.005 


±0.012 


0.920 ± 0.386 


0.43 


±0.36 


0.552 


± 


0.131 


2 


± 


17 


0.661 


± 


1.211 


2v 






156 


± 


71 


0.950 


± 0.072 


0.078 


±0.012 


0.996 ± 0.026 


0.93 


± 0.35 


0.391 


± 


0.060 


63 


± 


16 


0.956 


± 


0.042 


2v 


+ 


Viol 


108 


± 


74 


0.192 ±0.109 


0.006 


±0.012 


0.265 ± 0.302 


0.52 


±0.36 


0.067 


± 


0.108 


4 


± 


17 


0.235 


± 


0.598 


2v 


+ 


2v ro , 


193 


± 


74 


0.030 


± 0.061 


0.044 


±0.012 


0.023 ± 0.044 


1.50 


±0.36 


0.774 


± 


0.038 


35 


± 


17 


0.059 


± 


0.077 



The equivalent width W,\, radial velocity (V), the second 
moment (V 2 ) (linked to the line width), and the third moment 
(V 3 ) (a measure of the line asymmetry) were determined by 
integrating normalized spectra in the fixed 2 A window cen- 
tred at the Ndm line. The relative precision of the RV mea- 
surements is considerably higher compared to and other 
moments. Consequently, we have chosen (V) for detailed fre- 
quency analysis. The 1 860 RV determinations for the Nd in line 
are illustrated in Fig.|4] The data are phased with the rotational 
ephemeris 

HJD = 2448312.23606 + Ex 2.851976 (1) 

(Kurtz et al. U997t which are used throughout our paper. An in- 
homogeneous Nd distribution causes a prominent double-wave 



rotational modulation of RV. Pulsational variability is superim- 
posed on these longer term changes and manifests itself as a 
low amplitude high frequency scatter of the points in Fig.@] 

In the frequency domain (Fig.|5t) a contribution of non- 
radial pulsation is clearly seen as an excess of amplitude in the 
1350-1500 /jHz frequency range. In this paper we are primar- 
ily concerned with the pulsational changes of the spectrum of 
HR3831 and, hence, remove rotational modulation of RV by 
subtracting a 5-order Fourier fit (the smooth curve in Fig.@} 
from the data in time domain. The amplitude spectrum af- 
ter this prewhitening procedure is shown in Fig.|5j5. The high- 
amplitude peaks at low frequencies are gone and the pulsational 
contribution to the RV variability stands out with an unprece- 
dented clarity. Substantial pulsational signature in the (V) am- 
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Fig.5. Frequency analysis of the RV variation of the Ndm 6145.07 A line in the spectrum of HR3831. The left panels show 
low -resolution amplitude spectra. The vertical dashed lines mark the fundamental pulsation frequency (v = 1427.9958 yuHz) and 
its first harmonic (2v = 2855.916^Hz). a) The amplitude spectrum computed from the raw radial velocities shown in Fig.0J 
b) the amplitude spectrum after prewhitening of the rotational modulation. The right panels present high-resolution amplitude 
spectra in the vicinity of the main pulsation frequency. The vertical dashed lines show the expected position of the triplet around 
the main frequency: v— v rot , v, v+ v rot (where v rot = 4058.265 nHz). c) The amplitude spectrum after prewhitening of the rotational 
modulation; d) the corresponding CLEANed amplitude spectrum. 



plitude spectrum is observed only close to the main oscillation 
frequency of HR 3831. This differs from the photometric pulsa- 
tional behaviour (Kurtz et al. 1997 1 where significant amplitude 
excess is also seen at the first and second harmonics. 

The high-resolution amplitude spectrum (Fig.|5j;) shows 
that aliasing remains a problem in our single-site data de- 
spite a superb quality of the observational material. In a sit- 
uation where direct identification of the pulsation frequen- 
cies is hampered by the unfavourable window function, the 
CLEAN technique (Roberts et al. 119871 often aids in distin- 
guishing intrinsic pulsation frequencies from their aliases. We 
applied the CLEAN algorithm to the (V) measurements using 
200 iterations with gain=0.6. The resulting amplitude spectrum 
(Fig.|5jl) shows two high-amplitude peaks and a few weaker 
ones. Modification of the CLEANing parameters affects the 
low-amplitude peaks but the two main frequencies remain un- 
changed and, thus, can be identified without any ambiguity. 
Their position (respectively, 1424.06 and 1432.03 fjHz) coin- 
cides with the v - v rot and v + v lot variation seen in the high- 
speed photometry. This observation provides a strong argument 



for the overall similarity of the frequency structure of the pul- 
sational light and RV variation of HR 3831. 

Given the satisfactory agreement between the main features 
of the spectroscopic and photometric pulsational variability at 
the fundamental frequency of HR3831, we proceed with the 
least-squares analysis assuming that the oblique rotator model 
is valid for describing the (V) variation and that frequencies 
are always separated by multiples of v rot from the fundamen- 
tal frequency v and its harmonics. We use a non-linear least- 
squares algorithm to fit the whole fundamental septuplet and 
the first-harmonic quintuplet. (It should be emphasized that the 
respective 12 frequencies do not belong to physically differ- 
ent non-radial modes but arise from the rotational modulation 
of the same distorted oblique pulsation.) The radial velocity is 
approximated with the expression 

<V> = E A; cos {In ((r - 7o)(v + i v rot ) + ^)]+ 

' = 2 3 (2) 
£ 2 Aj cos [2n ((r - 7o)(2v + j v rot ) + <pj)], 

where the epoch of phase zero is To = 2448312.23606 and all 
amplitudes A (J , phases <fij, and main pulsation frequency v are 
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Fig. 6. Radial velocity curve for all observations of the Ndm 6145.07 A line in the spectrum of HR 3831. Symbols show individual 
RV measurements, with the rotational modulation of RV removed using the Fourier fit illustrated in Fig. [4] The solid curve shows 
RV predicted by the frequency solution discussed in the text. Individual panels correspond to the data obtained during each of the 
six observing nights. The number in the lower left of each panel gives HJD -2451900 of the start of spectroscopic monitoring. 
The x-axis is labeled with UT (in hours) and rotation phases. 



optimized simultaneously. The search for the best-fit rotation 
frequency can also be carried out. The resulting value of the 
stellar rotation period, P mt = 2.8624 + 0.0046 d, agrees with 
P rot = 2.851976 ± 0.00003 d determined by Kurtz et al. GTW7\ 
but is certainly less precise due to a shorter time-span of our ob- 
servations. Hence, we prefer to adopt the constant photometric 
Vrot = 4058.265 nHz in the fit with Eq. 0. The final parame- 
ters of the 12-frequency solution are reported in Table[2] Some 



of the pulsational contributions have amplitudes below the for- 
mal detection threshold but are, nevertheless, included here for 
completeness and to get an estimate of the upper limit of the 
RV variability at the respective frequencies. On the other hand, 
the second-harmonic triplet was not included in the fit because 
none of its components has amplitude above 16 ms _1 . 

The fundamental pulsation frequency determined in the 
least-squares analysis is v = 1427.9958 ± 0.0072 //Hz. This 
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Fig. 7. Amplitude spectra of the pulsational variability of (in this order, from top to bottom) the equivalent width Wa and the 
first three moments - the radial velocity (V), the second moment (V 2 ), and the third moment (V 3 ) - measured for the Ndm 
6145.07 A line. The left panels show low-resolution amplitude spectra. The vertical dashed lines marks the main pulsation fre- 
quency and its first harmonic. The right sequence of panels illustrate amplitudes of the components of the fundamental septuplet 
(v - 3v ro t, . . . , v, . . . , v + 3v rot ) derived with least-squares analysis. The horizontal dashed line corresponds to the 3<x uncertainty 
of the fitted amplitudes. The lower amplitude spectra in the left panels show the amplitude spectra after the frequency solutions 
illustrated in the respective right panels (see also Table[2) have been removed from the data. These residual amplitude spectra are 
shown at the same scale as the original spectra and are shifted downward for display purpose. 



is the best determination of the main pulsation frequency of 
HR 3831. The comparison between observed RVs of the Nd m 
line and prediction of the frequency solution is illustrated in 
Fig-IU Pulsational RV variation reaches 3.6 kms 1 peak-to- 
peak amplitude which is the highest value ever observed for 
a roAp star. For all observing nights the stellar pulsational 
behaviour is modelled with a good precision. The fundamen- 
tal septuplet plus the first-harmonic quintuplet explain 86% of 



the data variance. The residual amplitude spectrum of (V) is 
shown in the upper panel of Fig.0 The average noise level is 
14 ms -1 , whereas the highest high-frequency residual peaks 
reach ^40 ms~'. In the low-frequency domain we find a 
59 ms" 1 peak at v = 10.239 d -1 and its alias at v = 1 1.238 d" 1 . 
Detection of this variation is not statistically significant, but, if 
not caused by the instrumental effects, it can be a candidate for 
a weak 6 Scuti-type pulsation mode. 
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Fig. 8. Rotational modulation of the pulsation amplitude (left panels) and phase (right panels) of the variability of moments of 
the Ndm 6145.07 A line in the spectrum of HR3831. The panels show (in this order, from top to bottom) measurements of the 
equivalent width Wa, radial velocity (V), the second moment (V 2 ), and the third moment (V 3 ). Individual points in these diagrams 
were calculated by fitting the main pulsation frequency (v = 1427.9958 ;uHz) to the subsets of 38 spectroscopic observations 
(about 5 pulsation cycles or 56 min). The smooth curves show rotational modulation expected from the frequency solution for 
the fundamental septuplet given in Tableland illustrated in Fig.0 



Although we found a good overall concordance between 
the photometric and spectroscopic pulsation behaviour, the am- 
plitudes and phases of the individual frequency components 
turn out to be quite dissimilar. The first striking discrepancy 
between pulsational light variation and RVs is a significantly 
more sinusoidal behaviour of the latter observable. In fact, 
among the first-harmonic quintuplet components only the 2v 
variation is detected reliably and its amplitude is 14 times lower 
than the amplitude of the largest component of the fundamen- 



tal septuplet. Compare this with the time-resolved photometry 
in a Johnson B filter (Kurtz et al. 119971 for which this ratio is 
5 and even the second harmonic can be detected. Another dif- 
ference is a substantially weaker main frequency component 
observed in the RV variation. A widely used diagnostic of the 
oblique pulsator geometry (A + i + A_i)/Ao = tan i tan 6 (where 
i is inclination angle the rotation axis and S is the mode obliq- 
uity relative to this axis, see Kurtz & Martinez 2000 1 is 9.3 
according to Kurtz et al. ( 1997 1 but this parameter is 66.9 for 
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our spectroscopic data, indicating that one or both characteris- 
tic angles are close to 90°. The asymmetry between the ampli- 
tudes of the v + v rot pulsation components is also larger for RV: 
= 1.82 versus 1.21 in photometry. 

4.3. Moment variation and rotational modulation 

We proceed with analysis of the pulsational variation of 
the equivalent width and higher-order moments of the Nd in 
6145.07 A line using approach similar to the one applied to 
RV oscillations. The amplitudes and phases of the fundamental 
septuplet and first-harmonic quintuplet were optimized simul- 
taneously. The main pulsation frequency was kept constant at 
the value of v = 1427.9958 //Hz determined from the first mo- 
ment. Results of the fit of Wx, (V 2 ), and (V 3 ) are presented in 
Tableland are illustrated in Fig. [7] 

It is clear that pulsational variation of all line profile observ- 
ables considered here are detected with extremely high S /N at 
the main pulsation frequency. Weak harmonic variability is also 
found at ^4cr level for the (V 2 ) and (V 3 ) moments. Inspection 
of the residual amplitude spectra of the line profile moments 
(Fig. [7J does not reveal any additional statistically significant 
amplitude peaks. Therefore, the frequency model defined by 
Eq. (0 is fully adequate and explains all variability present in 
the time series of the line profile moments. 

An alternative presentation of the spectroscopic variability 
induced by the oblique non-radial pulsation in HR3831 can 
be obtained by considering pulsation with only the main fre- 
quency v but accounting for the rotational modulation by us- 
ing a time-dependent amplitude and phase. This formulation 
closely corresponds to the actual pulsation behaviour thought 
to be present in HR3831 and is mathematically equivalent to 
fitting constant amplitudes and phases of the fundamental sep- 
tuplet components. Rotational modulation of the pulsations in 
Wx, (V), <V 2 >, and (V 3 ) is illustrated in Fig. |8] The measure- 
ments shown in this figure were obtained by fitting cosine func- 
tion to 47 subsets of our observations covering roughly 5 pulsa- 
tion cycles. Rotational modulation of the line profile character- 
istics is well defined and is in perfect agreement with the modu- 
lation expected for the fundamental septuplet solution obtained 
above, confirming validity of both approaches to time-series 
analysis of the spectroscopic observations of HR 3831. 

The amplitude of Wx variation reaches w 5 mA at the ro- 
tation phase of the magnetic extrema. The phase modulation 
of the equivalent width changes is characterized by a sharper 
n radian transition at the crossover phase compared to the (V) 
phase modulation. Thus, Wx shows more similarity to the pho- 
tometric pulsational behaviour and, to some extent, can be used 
as a proxy of the latter. 

The radial velocity amplitude modulation is qualitatively 
similar to the Wx amplitude changes, but does not show asym- 
metry between the two pulsation poles evident in the Wx data. 
We note that amplitude changes of the RV variation found in 
our paper are consistent with the predictions of the generic 
oblique pulsator model. On the other hand, we do not confirm 
a single-wave rotational modulation of the REE RV amplitude 
claimed by Balona (2002i. Considering the large scatter of the 



measurements and a poor rotation phase coverage evident in 
Fig. 9 of Balona's paper, we conclude that no rotational modu- 
lation of the RV oscillations was detected in his study and the 
single-wave sinusoidal fit is not significant. 

The second moment amplitude shows an unexpectedly 
complex rotational modulation with the two maxima and two 
minima during each rotation cycle. The minima observed at the 
phases of maximum RV are deeper and narrower compared to 
the decrease of (V 2 ) amplitude associated with the pulsational 
equator (rotation phases ipx§25 and 0.75). Below we will 
demonstrate that this interesting behaviour of (V 2 ) is closely 
related to the distortion of oblique pulsation by the dipole com- 
ponent of the stellar magnetic field. 

The amplitude maxima in the rotational modulation curves 
of Wx and (V) presented in Fig.[8]occur slightly later than pre- 
dicted by the ephemeris of Kurtz et al. J1997i . This Aip = 
0.036 + 0.008 discrepancy between the phases of the photo- 
metric and spectroscopic pulsation maxima may be intrinsic to 
the star or indicate an error in the adopted rotation period. If 
the latter is true, P rot = 2.852055 ± 0.000017 d would bring 
in agreement the time of the photometric pulsation maximum 
(HJD=24483 12.23606, Kurtz et al.) and our spectroscopic ob- 
servations. The improved rotation period deviates by less than 
3<x from the P rot derived by Kurtz et al. All results presented 
in our paper are based on the data covering just 3.6 stellar rota- 
tions and, therefore, are unaffected by this small correction to 
the rotation period. 

5. Pulsational changes in other lines 

Observations of HR3831 analysed in this paper for the first 
time permitted measurements of pulsation amplitudes substan- 
tially below 50 ms -1 for individual spectral lines. Precision of 
our data enables detection of pusational RV variation for essen- 
tally all strong and medium-strength lines. Analyses of these 
absorption features was carried out using the same procedure 
as was outlined above for the Ndm 6145.07 A line. Table 
reports amplitudes of the components of the fundamental sep- 
tuplet for 6 lines and blends in the 6120-6150 A spectral inter- 
val. Especially important is the first definite detection and mea- 
surement of the pulsational RV variation for the light and iron- 
peak elements, represented in the studied spectral region by 
the Cai 6122.22 A and Fen 6147.74 A lines. The substantially 
larger amplitude of pulsations found for the Fen 6149.26 A 
line reflects blending by an unidentified REE spectral line (cf. 
Kochukhov & Ryabchikova 2001ai. Oscillations are also read- 
ily detected for the Ban 6141.71 A line, whose pulsation am- 
plitude in HR3831 is comparable to the pulsational changes of 
this line in the roAp star yEqu (Kochukhov & Ryabchikova 
I2()()lal . 

Our measurements of the RV variation of metal lines 
(Table and the Ndm 6145.07 A spectral feature (Table 
show a factor of 50 difference in amplitude and significant 
(sometimes as large as n radian) discrepancy in phase. This 
picture is consistent with general explanation of the pulsational 
spectroscopic behaviour of roAp stars proposed by Kochukhov 
& Ryabchikova (2001a) and Ryabchikova et al. (2002 1: the 
running pulsation wave propagates outward through the chemi- 
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Fig.9. Amplitude spectra of the RV variation observed for spectral features other than the Ndm 6145.07 A line. The vertical 
dashed lines mark the fundamental pulsation frequency (v = 1427.9958 /jHz) and its first harmonic (2v = 2855.916 //Hz). 



Table 3. Results of non-linear least-squares fit of the RV variation of spectral lines other than Nd in 6145.07 A. 



Frequency 


Cai 6122.22 A 


Cm 6129.22 A 


Sii + Cen 6131.9 A 




A (ms 


- 1 ) 


V 


A (ms- 1 ) 




f 


A (ms- 1 ) 




v - 3v rot 


25 : 


t9 


0.511+0.055 


44+19 


0.509 


± 0.068 


23 ±29 


0.484 ± 0.206 


v - 2v rot 


12: 


t8 


0.446 + 0.108 


24+18 


0.289 


±0.116 


45 ±27 


0.036 ± 0.098 


v - v IO , 


62: 


t8 


0.439 + 0.021 


67 ± 17 


0.211 


±0.041 


49 ±27 


0.470 ± 0.087 


V 


27 : 


t9 


0.095 + 0.051 


21 ± 18 


0.807 


±0.143 


35 ±29 


0.189 ±0.131 


V + Vrot 


53: 


t8 


0.157 + 0.024 


115 ± 17 


0.211 


± 0.024 


81 ±27 


0.117 ±0.053 


v + 2v rot 


10: 


t8 


0.407 + 0.133 


28 ± 18 


0.674 


±0.100 


38 ±27 


0.288 ±0.113 


v + 3v ro , 


13: 


t9 


0.108 + 0.107 


43 ± 19 


0.099 


± 0.069 


51 ±29 


0.206 ± 0.092 


Frequency 


Ban 6141.71 A 


Fen 6147.74 A 


Fen + REE 6149.26 A 


v - 3v rot 


41 ± 


11 


0.326 ± 0.042 


41 ± 10 


0.536 


± 0.040 


86 ± 12 


0.392 ± 0.021 


v - 2v rot 


14 + 


10 


0.239 + 0.122 


15 ± 10 


0.591 


±0.105 


47 ± 11 


0.144 ±0.037 


V ~ Vrot 


132 ± 


10 


0.302 + 0.012 


45+9 


0.573 


± 0.033 


248 ±11 


0.287 ± 0.007 


V 


16 ± 


11 


0.018 + 0.109 


19 + 10 


0.895 


± 0.085 


18 ± 11 


0.779 + 0.103 


V + Vrot 


81 ± 


10 


0.173 +0.020 


47+9 


0.093 


± 0.032 


146 + 11 


0.985 +0.012 


v + 2v rot 


26 ± 10 


0.567 ± 0.064 


4+10 


0.530 


±0.398 


60 ± 11 


0.563 ± 0.029 


v + 3v rot 


8 + 


11 


0.269 + 0.210 


25 ± 10 


0.205 


± 0.064 


61 ± 12 


0.007 ± 0.030 



cally inhomogeneous stellar atmosphere with increasing ampli- 
tude. Large difference in the pulsation characteristics of various 
lines arises from an interplay between the depth-dependence of 
magnetoacoustic wave and chemical gradients formed by ra- 
diative diffusion. In particular, the rare-earth elements are con- 
centrated in high atmospheric layers and probe regions with 
a large pulsation amplitude, whereas Ca and Fe lines origi- 
nate from the deeper layers where pulsations are very weak. 
The strong Ban lines are probably sensitive to the intermediate 
depths. 



6. Interpretation of moment variations 

The problem of modelling variation of the line profile moments 
for stars pulsating in oblique non-radial modes was addressed 
for the first time in a recent paper by Kochukhov (2005 1. In 
this study we derived relevant expressions and developed a nu- 
merical procedure suitable for analysis of the time-resolved 
spectroscopic observations of roAp stars. It was demonstrated 
that a non-axisymmetric superposition of the pulsation and ro- 
tation velocity fields in oblique pulsators leads to a qualita- 
tively new behaviour of some of the line profile characteristics. 
Furthermore, a distortion of the oblique dipole modes predicted 
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Fig. 10. Comparison between the observed rotational modulation of the pulsation amplitude and phase of the (V) and (V 2 ) 
moments of the Ndm 6145.07 A line (symbols) and predictions of the oblique pulsator model discussed in the text. The solid 
lines shows results obtained with the best-fit combination of the non-axisymmetric dipole (( = l,m = -1, 0, 1) and axisymmetric 
octupole {( = 3, m = 0) pulsation components, whereas the dashed line illustrates moment modulation expected from the dipole 
components alone. 



by the recent theoretical studies of the stellar magnetoacoustic 
oscillations (Saio & Gautschy 2004; Saio 2005) can be readily 
diagnosed through the moment analysis. In particular, the shape 
of the pulsation phase modulation for the radial velocity and 
second moment is very sensitive to non-axisymmetric pulsa- 
tion components, whereas the rotational modulation of the sec- 
ond moment amplitude is best suited to revealing axisymmet- 
ric magnetically induced distortion of pulsations. The present 
paper deals with the first application of the new moment tech- 
nique to observations of a real roAp star and aims at deriving a 
detailed model of the stellar pulsation velocity field. 

The moment code of Kochukhov ( 2005 ) was employed to 
compute time-dependent amplitudes and phases of (V) and 
(V 2 ) for 100 equidistant rotation phases. The local line pro- 
files were approximated by a Gaussian with <x = 5 kms~'. 
Inclination angle i = 68.2° and projected rotation velocity 
v e sini = 33.0 kms~' were adopted following results of the 
model atmosphere and surface structure analysis of HR3831 
presented by Kochukhov et al. (|5t)04). The free parameters of 
our spectroscopic oblique pulsator model included the obliq- 
uity 5 and phase angle \ giving orientation of the pulsation axis 
and amplitudes of the I = 1, m = 0, ±1 and £ = 3, m = pul- 
sation components. These parameters were optimized through 
the comparison of the model predictions with the observed ro- 
tational modulation of the (V) and (V 2 ) amplitudes and the first 



moment phase. An equal weight was given to all three observ- 
ables and horizontal pulsation motions were neglected. 

The final set of parameters of the oblique pulsator model 
required to fit variation of the Ndm 6145.07 A line profile mo- 
ments is summarized in Tabled Comparison between observa- 
tions and model predictions is illustrated in Fig.^| All consid- 
ered pulsation observables except the phase of (V 2 ) are either 
reproduced successfully ((V) amplitude and phase) or fairly 
well ((V 2 ) amplitude) by our model. It appears that complex 
modulation of the second moment amplitude is explained by 
the presence of I — 3, m = pulsation. It definitely cannot be 
reproduced by any combination of the dipole pulsation compo- 
nents (see Fig.llQI alone. We find that at the formation depth 
of the Nd in 6145.07 A line pulsation geometry is slightly non- 
axisymmetric, with the amplitude of the £ = 1, m = ±1 compo- 
nents not exceeding ~ 15% of the contribution due to the £ = 1, 
m — spherical harmonic. On the other hand, we uncover a 
prominent distortion of the dipole pulsation in HR 3831 by the 
global magnetic field. It manifests itself through the presence 
of a substantial axisymmetric octupole {( = 3, m = 0) pulsa- 
tion component and leads to enhancement of velocity ampli- 
tude close to the pulsation poles. These findings are in good 
agreement with the outcome of pulsation Doppler imaging of 
HR3831 (Kochukhov 2004b) and are in line with theoretical 
predictions of Saio & Gautschy (2004i and Saio (2005). 
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Table 4. Parameters of the oblique pulsator model of HR 3831 
determined from the variation of the Nd m 6145.07 A line pro- 
file moments. This table reports stellar inclination i, angles 6 
and x specifying orientation of the pulsation axis, and ampli- 
tudes (renormalized to yield the maximum surface velocity) of 
different spherical harmonic contributions to the pulsation ve- 
locity field. 



Parameter Value 







i(l 


68.2 






SC) 


88.5 ±0.6 






xO 


167.1+0.5 


A(i = 


1 


,m = 0) (kms- 1 ) 


3.33 + 0.03 


A{{ = 


1, 


m = ±1) (kms -1 ) 


0.51+0.02 


A(l = 


3 


,m = 0) (kms- 1 ) 


2.90 + 0.11 



Precise information about orientation of the pulsation axis 
obtained in our study demonstrates that oblique oscillations in 
HR3831 are aligned with the dipolar magnetic field compo- 
nent. The latest and most detailed modelling of the magnetic 
topology of HR 3831 (Kochukhov et al. 2004 (has yielded mag- 
netic obliquity (3 = 86.8 + 6.2° which coincides with the pul- 
sation mode obliquity 6 = 88.5 + 0.6° determined through the 
moment analysis. 

The present model describing variation of the Ndm 
6145.07 A line cannot be improved by introducing high-order 
non-axisymmetric pulsation contributions. A formal least- 
squares analysis also yields the ratio of the horizontal to 
vertical pulsation amplitude, K, indistinguishable from zero. 
However, the latter result may reflect a failure of the standard 
constant-/^ parameterization of the horizontal pulsation mo- 
tions adopted in the moment technique. In fact, this approxima- 
tion is poorly suited to characterize latitude-dependent horizon- 
tal pulsation amplitude expected for magnetoacoustic pulsators 
(Saio & Gautschv l27704T i. 

The most significant limitation of the modelling described 
in this Section is related to the fact that we disregard in- 
homogeneous distribution of Nd at the surface of HR3831. 
Neodymium surface map reconstructed by Kochukhov et al. 
(2004;) shows relative underabundance of this element along 
the magnetic equator. Consequently, the moment diagnostic 
is somewhat less sensitive to the pulsation wave properties in 
these regions of the stellar surface. However, we do not expect 
results of the overall moment analysis to be significantly dis- 
torted by the presence of Nd spots because in our modelling 
information is mostly extracted from the high-amplitude LPV 
observed at the rotation phases when no large Nd abundance 
gradients are present at the visible stellar disk. 

7. Discussion and conclusions 

Despite a long history of the high-speed photometric obser- 
vations of pulsations in roAp stars, mechanisms responsible 
for the pulsational luminosity variation are not understood. In 
particular, no theoretical model has been able to explain the 
observed complexity (non-linearity, peculiar wavelength de- 
pendence, etc.) of the photometric time series. Thus, no sim- 
ple relation can be established between changes of luminos- 



ity and fundamental pulsation characteristics, such as the pul- 
sation displacement vector. This is why all interpretations of 
the photometric pulsational variability of roAp stars start with 
the assumption that geometrical picture of the relative lumi- 
nosity changes, 61 1 1, is the same as for the vertical component 
of the displacement vector and is given by a sum over spher- 
ical harmonic functions. Thus, even the latest state-of-the-art 
theoretical investigations into pulsational light variability (e.g. 
Bigot & Dziembowski 2002ji have to abandon hopes to build 
a physically consistent picture and present essentially a geo- 
metrical model of the continuum intensity variation associated 
with p-mode oscillations. A progress in this field requires de- 
tailed non-adiabatic MHD modelling (first models of this kind 
are being developed by Saio (2005l) as a basis for estimating 
contributions of different terms to the observed luminosity vari- 
ation of roAp stars. Nevertheless, the problem of a low intrinsic 
diagnostic content of a disk-integrated photometric observable 
(cf. Saio & Gautschy 2004) could never be overcome. 

In contrast to unsatisfactory understanding and limited 
information content of the pulsational light curves, spectro- 
scopic observations of rapid oscillations permit a straightfor- 
ward astrophysical interpretation. The observed pulsation ve- 
locity fluctuation is directly related to the time depedence of 
the pulsation displacement vector. Given the absence of signif- 
icant pulsational changes of the strength or equivalent width 
of diagnostic lines, we can also assert that the observed RV 
variation is not visibly influenced by secondary effects, such 
as rapid changes of temperature and density of atmospheric 
plasma. Consequently, interpretation of the spectroscopic ob- 
servations of roAp pulsations offers us a real opportunity to 
develop physical models of magnetoacoustic pulsation for the 
first time since the discovery of roAp stars. Furthermore, the 
presence of strong vertical gradients of chemical composition 
in the atmospheres of roAp stars enables an unprecedented ver- 
tical resolution of pulsation waves. This makes spectroscopic 
studies of roAp pulsations a very promising research direction, 
both in terms of its enourmous diagnostic potential and possi- 
bility of establishing a physically realistic model of the pulsa- 
tional phenomena from time-resolved observations. 

The line profile analysis has become a standard tool for the 
investigation into pulsational variability of many types of oscil- 
lating stars. Yet, due to their short pulsation periods, roAp stars 
were not studied before with this technique. The present report 
on LPV in the prototype roAp star HR3831 fills this gap and 
represents a significant step towards developing observational 
methods and analytical tools able to bring us closer to under- 
standing intricate nature of the magnetoacoustic pulsations in 
magnetic A stars. 

Our line profile study of HR3831 focuses on the vari- 
ability of the Ndm 6145.07 A line, which shows maximum 
amplitude of the pulsational variability among spectral fea- 
tures in the studied wavelength region. It also exhibits one of 
the largest pulsation amplitudes ever observed in a roAp star. 
Spectroscopic material described in our paper shows unam- 
bigous evidence of the pulsational line profile variability and 
presents the very first comprehensive study of this phenomenon 
in a star pulsating in oblique non-radial mode. In particular, it 
enabled us to detect and characterize rotational modulation of 
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Fig. 11. Depth dependence of the magnetoacoustic wave properties in the atmosphere of HR 3831. a) (A_i +A+0/2 characterizes 
variation with depth of the average amplitude of the v + v lot components, d) depth dependence of the respective phases of the 
v + v ro t (open circules) and v - v, ot (filled circles) pulsation, c) the (A_i + A+i)/Ao parameter characterizes strength of the main 
frequency component relative to the v+v rot components, d) A_i/A+i measures asymmetry of the v±v rot components. The symbols 
show results obtained for individual spectral lines, whereas the dashed line corresponds to the photometric pulsation behaviour 
in a Johnson B filter. 



the LPV pattern - an effect expected for roAp stars long ago 
(e.g. Baade & Weiss [1987) but never observed until now. 

Pulsational changes of the equivalent width and the first 
three moments of the Nd in line are detected with a very high 
S IN ratio and studied over the whole stellar rotation cycle. 
Unlike strongly non-sinusoidal pulsation light curves (Kurtz et 
al.QP.97l, all spectroscopic observables show only a marginal 
presence of the first and no signature of the second harmonic in 
their amplitude spectra. Following Balona (2002i, we interpret 
this observation to be an evidence that spectroscopy probes si- 
nusoidal pulsation of the star, whereas variation of temperature 
and other parameters contributing to luminosity fluctuation is 
related to the pulsation displacement in a complex non-linear 
way. 

Time-series analysis of the RV variation resulted in a re- 
markably accurate determination of the main pulsation fre- 
quency of HR3831 and possible improvement of the rotation 
period of the star. 

Taking advantage of the comprehensive observation of the 
LPV in the Nd in 6145.07 A line, we develop a physically con- 
sistent oblique pulsator model of HR 3831. This model directly 
characterizes horizontal geometry and amplitude of the pulsa- 



tion displacement in the upper atmospheric layers where strong 
REE lines are formed. Rapid variation of the first two line pro- 
file moments is interpreted using a new version of the mo- 
ment mode identification method (Kochukhov 2005 1. We de- 
rive orientation of the pulsation axis and amplitudes of vari- 
ous spherical harmonic contributions to the pulsation geometry. 
These results generally agree with the pulsation Doppler imag- 
ing map of velocity field obtained for HR3831 in our earlier 
study (Kochukhov 200415). 

Using the first detailed map of pulsations in a roAp star 
we are able to test predictions of recent theoretical studies of 
magnetoacoustic pulsations. Two main types of the distortion 
of the basic oblique dipole pulsator geometry were suggested. 
On one hand, Bigot & Dziembowski ( 2002 1 take into account 
the interaction between weak magnetic field, pulsation, and ro- 
tation and predict a non-axisymmetric dipole pulsation geom- 
etry, which is not aligned with the magnetic field axis. On the 
other hand, Saio & Gautschy (I2004t and Saio J2005> neglect 
the stellar rotation and instead suggest that the dominant dis- 
tortion of pulsation appears due to the strong magnetic field 
influence. Such a distortion manifests itself in the presence of 
substantial o&d-C harmonic components in the pulsation ge- 
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ometry. Our moment analysis results show a small (but sig- 
nificant) non-axisymmetric contribution, which could be as- 
cribed to the I = 1, m = +1 components, but we find no 
misalignment of the pulsation and magnetic field axes sug- 
gested by Bigot & Dziembowski ( I2UU21 for HR3831. At the 
same time, we reveal a strong I — 3, m — pulsation com- 
ponent, in a good agreement with the model developed by Saio 
& Gautschy ( 2004). This demonstrates that magnetic distortion 
of pulsation modes is by far more important effect compared to 
the rotation and magnetic field coupling investigated by Bigot 
& Dziembowski (2002). Although a future ultimate model of 
roAp pulsations should certainly include effects scrutinized by 
Bigot & Dziembowski, neglect of the most important magnetic 
perturbation of p-modes in their modified oblique pulsator the- 
ory makes it inapplicable to roAp stars with the average mag- 
netic field strength above su 1 kG, which includes HR 3831. On 
the other hand, success of the theoretical framework proposed 
by Saio & Gautschy (2004) demonstrates that their model is 
based on a realistic description of the physics of magnetoacous- 
tic pulsation and, hence, represents a major progress in under- 
standing excitation and propagation of pulsation waves in roAp 
stars. 

Unprecedented precision of the time-series spectroscopic 
observations of HR 3831 analysed here permitted the first def- 
inite detection of weak pulsational RV changes in the lines 
of light and iron-peak elements. This makes it possible to at- 
tempt a tomographic mapping (Kochukhov l2003l ) of the verti- 
cal structure of the pulsation mode in HR3831. Such a proce- 
dure requires accurate information about formation depths of 
the spectral lines involved in the analysis. Chemical stratifica- 
tion significantly affects the depths from which lines of indi- 
vidual ions originate and this is why a model atmosphere anal- 
ysis and stratification diagnostic are necessary prerequisites to 
reconstructing a realistic picture of the depth dependence of 
pulsation characteristics. Such a study lies outside the scope of 
this paper and cannot be performed using spectra with a limited 
wavelength coverage available to us. Nevertheless, it is possible 
to probe certain general properties of the pulsation wave struc- 
ture by using results of the tomographic mapping of pulsations 
in y Equ (Ryabchikova et al. 2002 ), which was largely based on 
the same set of spectral lines as analysed here. Similarity of the 
atmospheric parameters and chemical composition of HR 3831 
and y Equ suggests that the vertical abundance profiles and line 
formation depths are comparable in the two stars. Thus, we can 
assert that the Ndm 6145.07 A line is formed in a high atmo- 
spheric layer, the Ban 6141.71 A line is sensitive to interme- 
diate depths, whereas lines of Ca and Fe should originate from 
deeper layers. At the same time, the line of Cai 6122.22 A 
investigated here is considerably stronger that Fen 6147.71 A 
and, hence, it should form somewhat higher compared to the 
iron line. Using this tentative assignment of the relative for- 
mation depths of the four lines we are able to obtain the first 
picture of the vertical cross-section of the pulsation mode in 
HR 3831. Information about depth-dependent behaviour of the 
amplitudes and phases of the fundamental triplet is summarized 
in Fig.Gj] 

Our observations establish that the amplitude of pulsations 
increases by a factor of 17 from the lower to upper layers in the 



atmosphere of HR3831 (Figll lfa). Such a behaviour is an ex- 
pected consequence of the constancy (or only a small change) 
of the pulsation energy density pv 2 /2 in the presence of a rapid 
decrease of the gas density in stellar astmosphere. Large pul- 
sation phase difference between the RV variations of spectral 
lines sensitive to upper and lower atmospheric layers clearly 
points to a short vertical length of pulsation waves. Comparison 
of the pulsation phases derived for different lines in the studied 
spectral region of HR3831 does not reveal a bimodal phase 
distribution but shows a quasi-continuous phase change with 
depth (Figll lb). Therefore, we conclude that the vertical struc- 
ture of pulsation in the atmosphere of HR3831 resembles a 
running wave rather than a standing wave. With the data avail- 
able here, we find no evidence for the existence of pulsation 
node in the atmosphere of HR 3831. 

Comparison of our results with the theoretical non- 
adiabatic calculation (Saio 2005) shows a reasonable agree- 
ment with respect to the magnitude of the pulsation amplitude 
increase towards the upper atmospheric layers and the absence 
of pulsation node in the line-forming region. However, a signif- 
icant depth dependence of pulsation phase is not expected for 
this theoretical model. 

In addition to probing depth dependence of the average pul- 
sation amplitude and phase, our measurements of the RV vari- 
ation of the lines of different ions enable extraction of the in- 
formation about changes of the pulsation wave geometry with 
height. A common procedure to characterize oblique pulsator 
geometry and its distortion by the effects of stellar rotation is 
to study asymmetry of the v + v rot peaks and estimate relative 
amplitude of the variation at the unshifted frequency v. The 
respective observables, A-i/A+\ and (A_i + A + i)/Aq, derived 
from the light curves of roAp stars have been interpreted in 
a number of theoretical investigations (see Kurtz & Martinez 
2000 for recent overview) under the assumption that these pa- 
rameters permit a simple geometrical interpretation within the 
oblique pulsator framework. Our study of HR 3831 reveals that, 
similar to a dramatic depth dependence of the pulsation ampli- 
tude and phase, relative amplitudes of the triplet components 
do not have well-defined values consistent with those inferred 
from the photometric observations of HR 3831 but change sig- 
nificantly with height in the stellar atmosphere. In particular, 
analysis of the Fe n and Ca i lines shows that pulsation in the 
lower part of the photosphere has a moderate asymmetry of the 
v + Vrot components and noticeable contribution of the pulsa- 
tion at the main frequency v. Both characteristics agree with 
photometric results (Kurtz et al. 1997). However, asymmetry 
of the rotationally separated peaks increases and pulsation at 
the unshifted frequency decreases below the formal detection 
threshold when the diagnostic Nd m line sensitive to the upper 
atmospheric layers is considered (see Fig.ll It. d). 

Observed pulsation behaviour of different lines in the 
HR3831 spectrum can be attributed to a signature of the depth 
dependent superposition of the acoustic and magnetic pulsa- 
tion wave components. Importance of the latter contribution 
presumably increases with height following the correspond- 
ing increase of the ratio of magnetic to gas energy density. 
Furthermore, in the lower part of atmosphere dissipation of the 
magnetic slow waves and related horizontal pulsation motions 
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(Saio & Gautschy 2004 1 could also contribute to the observed 
spectroscopic and photometric pulsation variability, introduc- 
ing further difference with oscillation observed in the upper at- 
mosphere using REE lines. Thus, a common notion of A_i/A + j 
and (A_i + A +1 )/Aq as "geometrical" parameters typical of the 
whole pulsation structure is meaningless. Both quantities show 
significant variation with depth and should always be associ- 
ated with specific layers of the atmosphere of a pulsating mag- 
netic star. Thus, extreme caution is called for a purely geomet- 
rical interpretation of the relative amplitudes of the v,v ± v ro t 
components. 

Interpretation of the unique high time resolution spectro- 
scopic data presented here inevitably calls for a significantly 
more complex three-dimensional picture of the structure of pul- 
sating cavity in roAp stars compared to the plane-parallel view 
provided by the traditional oblique pulsator paradigm. There 
is clearly a need to advance from a rudimentary phenomeno- 
logical model to a self-consistent physical description of the 
propagation of magnetoacoustic waves in stellar atmospheres. 
The first important step in this direction is self-consistent mod- 
elling of the vertical and horizontal chemical inhomogeneities 
and pulsational variability using pulsation Doppler Imaging 
method (Kochukhov 2004a 2004b I. This modelling can be ap- 
plied to spectral lines formed at different optical depths in order 
to deduce a 3-D structure of pulsation velocity field. This mod- 
elling is currently underway for HR 3831 and will be presented 
in a forthcoming publication. 
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